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ABSTRACT 

In this work, we present results for the photometric and clustering properties of 
galaxies that arise in a A cold dark matter hydrodynamical simulation of the local uni- 
verse. The present-day distribution of matter was constructed to match the observed 
large scale pattern of the IRAS 1.2-Jy galaxy survey. Our simulation follows the for- 
mation and evolution of galaxies in a cosmological sphere with a volume of ~ 130 3 
Mpc 3 including supernova feedback, galactic winds, photoheating due to an uniform 
meta-galactic background and chemical enrichment of the gas and stellar populations. 
However, we do not consider AGNs. In the simulation, a total of ~ 20000 galaxies 
are formed above the resolution limit, and around 60 haloes are more massive than 
~ 10 14 M Q . Luminosities of the galaxies are calculated based on a stellar population 
synthesis model including the attenuation by dust, which is calculated from the cold 
gas left within the simulated galaxies. Environmental effects like colour bi-modality 
and differential clustering power of the hydrodynamical galaxies are qualitatively sim- 
ilar to observed trends. Nevertheless, the overcooling present in the simulations lead 
to too blue and overluminous brightest cluster galaxies (BCGs). To overcome this, we 
mimic the late-time suppression of star formation in massive halos by ignoring recently 
formed stars with the aid of a simple post-processing recipe. In this way we find lumi- 
nosity functions, both for field and group/cluster galaxies, in better agreement with 
observations. Specifically, the BCGs then follow the observed luminosity-halo mass 
relation. However, in such a case, the colour bi-modality is basically lost, pointing 
towards a more complex interplay of late suppression of star formation than what is 
given by the simple scheme adopted. 

Key words: hydrodynamics - galaxies: formation - evolution - cosmology: theory - 
methods: numerical 



1 INTRODUCTION 

The observational study of galaxy populations has seen an 
outstanding progress in recent years. With the advent of 
large galaxy redshift surveys, such as the Sloan Digital Sky 
Survey (SDSS; York et al. 2000) and the Two-degree Field 
Galaxy Redshift Survey (2dFGRS; Colless et al. 2001) it was 
possible to extend our knowledge of the local Universe to a 
new level of accuracy. 

In particular, it has been possible to carry out a robust 
determination of the luminosity function (LF) for galaxies 
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in the field (Norberg et al. 2002; Blanton et al. 2003) in 
different spectral bands and to better establish it for those 
galaxies populating denser environments, such as groups and 
clusters (e.g. Popesso et al. 2004, 2006). Consistently with 
previous work (e.g. Lin et al. 1996; Colless et al. 1999) it 
has been found that the field LF is well described by a sin- 
gle Schechter function requiring, however, a higher value for 
the luminosity density of the Universe indicating that pre- 
vious surveys suffered from selection effects and systematics 
due to photometry (Blanton et al. 2001). For higher density 
environments, and by means of the RASS-SDSS galaxy clus- 
ter survey, Popesso et al. (2004, 2006) found that a double 
Schechter function better fits the composite cluster LF ac- 
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Figure 1. Full-sky map of the simulated local volume in supcrgalactic coordinates. The map displays the projected gas density distribution 
up to a distance of ~ 80 h" 1 Mpc from the observer where main galaxy clusters are shown. Note that Virgo cluster, being the closest 
one to the observer, is particularly prominent. 



counting for a possible upturn in the number of faint objects 
as the luminosity of member galaxies decrease. 

Moreover, the large number of galaxies with measured 
spectroscopic redshifts, compared to that obtained in the 
past, has made possible the determination of their two-point 
clustering properties in a reliable way out to scales of order 
of tens of Mpc. These studies have shown that the real, pro- 
jected and redshift-space correlations can be well described 
by a decreasing power-law function that depends on the sam- 
ple colour, having a correlation length that increases with 
absolute magnitude (e.g. Norberg et al. 2001, 2002; Zehavi 
et al. 2002; Madgwick et al. 2003; Hawkins et al. 2003). 

Both the LFs and the clustering properties of the galax- 
ies are tools of fundamental importance in the study of 
galaxy formation since they provide a way to describe the 
most basic galaxy statistics. A successful model for structure 
formation must account for these observations. 

From the theoretical point of view, the widely spread 
semi- analytic models (SAMs) of galaxy formation (e.g. 
Kauffmann et al. 1999; Springel et al. 2001; Mathis et al. 
2002; De Lucia, Kauffmann & White 2004; Springel et al. 
2005; Bower et al. 2006; Cattaneo et al. 2006; Croton et 
al. 2006; Lagos, Cora & Padilla 2008; Fontanot et al. 2009; 
Guo & White 2009) provide a way to study the properties of 
galaxy populations with the advantage of a rapid exploration 
of the parameter space. In this approach, the galaxy popu- 
lation is followed within the skeleton provided by a parent 
dark matter simulation with the aid of physically motivated 
recipes to describe the different baryonic processes involved. 
These studies have pointed out the need of limiting excessive 
gas condensation in massive haloes to avoid the formation of 
very bright central galaxies (also known as brightest cluster 
galaxies; BCGs) that are inconsistent with observation. The 



usually invoked channel responsible for the star formation 
(SF) quenching in massive haloes is the active galactic nu- 
cleus (AGN) phenomenon. Similarly, it has been noted that 
in order to prevent an excess luminosity in the faint end of 
the galaxy field LF, a relatively strong SN feedback would 
be needed in smaller systems. 

On the other hand, cosmological hydrodynamical sim- 
ulations of galaxy samples, without including the effects of 
AGN feedback, have also been used to study the building 
up of the structure and their resulting properties in periodic 
boxes (e.g. Pearce et al. 2001; White et al. 2001; Yoshikawa 
et al. 2001; Oppenheimer & Dave 2006, 2008; Dave & Op- 
penheimer 2007; Ocvirk et al. 2008). In a similar way, using 
the dubbed zooming technique, several authors simulated 
hydrodynamical galaxies within high-resolution regions in a 
cosmological framework (e.g. Dolag et al. 2005; Saro et al. 
2006, 2008, 2009; Crain et al. 2009). In particular, Saro et 
al. (2006) resimulated a set of galaxy clusters studying the 
composite cluster LF and the environmental induced galaxy 
properties in these high density regions. These authors show 
that it is possible to reproduce the general observed trends 
of the cluster galaxy population, including, e.g. the colour- 
magnitude relation, the age and colour cluster-centric dis- 
tance dependence and the cluster LF (although the bright- 
end is affected by the presence of very bright BCGs). 

The present work makes use of the same zooming tech- 
nique to simulate the observed local volume following the 
chemical enrichment of gas and stars until the present epoch. 
This enables us to estimate in a consistent way the lumi- 
nosities and colours of galaxies formed in hydrodynamical 
simulations. The simulated volume is big enough to com- 
pute LFs, both for field and group/cluster galaxies, as well 
as galaxy correlation functions. In order to infer the way 
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Figure 2. Estimated optical depth profile in the V-band for one 
of the most massive BCGs in the simulation including extended 
intra-cluster light (M* ~ 10 13 Mq) averaged over radial pro- 
jected bins (top panel) and its correspondingly smoothed attenu- 
ation fraction due to dust (bottom panel) for the bj- and isf-bands 
(open and filled circles respectively). 

feedback acts in massive haloes (e.g. AGN feedback), we 
use a similar approach as SAMs, applying a simple recipe 
to the post-processed data of the simulation at z = 0. In 
this way, we ignore a fraction of the late-formed stars which 
would not appear in such a scenario and study the resulting 
effects on the luminosity-dependent properties of the simu- 
lated galaxies. It is important to note that, within the simple 
post-processing scheme adopted here, the simulations are no 
longer fully self-consistent, as also stars which are quenched 
by our post-processing procedure still interact with the sur- 
rounding medium. However, this effect is not expected to 
change our results significantly and can only be accessed 
with the next generation of hydrodynamical simulations in- 
cluding sub-scale models for AGN feedback more directly. 

The paper is organized as follows. In Section [2] wc 
describe the hydrodynamical cosmological simulation, to- 
gether with the method used to compute galaxy luminosi- 
ties and the associated dust-obscuration. In Section [3] we 
show the main results, presenting luminosities (for field and 
group/cluster galaxies), galaxy colours and correlation func- 
tions. We also discuss the implemented recipe to suppress SF 
in massive haloes. Finally, we close the paper with a sum- 
mary and our conclusions in Section [4] 



2 THE SIMULATION 

The cosmological simulation analyzed in this paper was gen- 
erated to reproduce the large scale distribution of matter in 
the local universe within the context of a flat ACDM sce- 
nario. The cosmological parameters at present time are a 
matter density parameter fi m = 0.3, a baryon density pa- 



Table 1. Simple SF suppression recipe: z cu t = (V max — Vq)/Vi, 
where z cu t represents the redshift since which we do not consider 
star formation in a given galaxy. Systems having V ma x < Vo are 
not affected by the suppression. 



Model 


Vo 


Vl 




[km s _1 ] 


[km s" 1 ] 
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100 


150 
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100 


200 


3 


100 


250 
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150 


150 


5 


150 


200 


6 


150 


250 


7 


200 


150 


8 


200 


200 


9 


200 


250 



rameter Sib = 0.04, a Hubble constant Ho = 100 h km s" 1 
Mpc with h = 0.7, and a rms density fluctuation erg = 0.9, 
which correspond to the Wilkinson Microwave Anisotropy 
Probe (WMAP) 1 year best-fitting cosmology (Spergel et 
al. 2003). 

In order to reproduce the local matter density field, the 
galaxy distribution in the IRAS 1.2-Jy galaxy survey (Fisher 
et al. 1994, 1995) was Gaussianly smoothed on a scale of 7 
Mpc and then linearly evolved back in time up to z — 50 
with the method proposed by Kolatt et al. (1996). The re- 
sulting high redshift field is then used as a Gaussian con- 
straint (Hoffman & Ribak 1991) to assign the unperturbed 
positions of 2 x 51 million gas and dark matter particles 
that are arranged on a glass-like distribution (e.g. Baugh, 
Gaztanaga & Efstathiou 1995). The mass of gas and dark 
matter particles is set to 4.7 x 10 8 and 3.1 x 10 9 h" 1 M 
respectively. Every gas particle is able to produce three gen- 
erations of stars, which results in a star particle mass of 
~ 1.6 x 10 8 /i" 1 M . 

We follow the evolution of the collisionless and gas com- 
ponents within a sphere with a comoving diameter of 160 
h^ 1 Mpc that is embedded in a periodic cosmological box 
of 340 Mpc on a side. This makes the volume of the 
simulated universe a factor of ~ 2 times larger than in pre- 
vious studies (e.g. Pearce et al. 2001). The region outside the 
sphere is sampled with ~ 7 million low resolution dark mat- 
ter particles, that allow us to asses the effect of long-range 
gravitational tidal forces. The comoving softening length for 
gravitational forces is set to 7 h^ 1 kpc (Plummer equiva- 
lent), value that is similar to the average particle separation 
in the cores of the densest simulated clusters. In Fig. [T] we 
show the projected gas density distribution of the simulated 
local sphere. 

The simulation was run using GADGET-2 (Springel 
& Hernquist 2002, 2003; Springel 2005) which is a Tree- 
smoothed particle hydrodynamics (SPH) code that fully 
conservates entropy during the evolution of the gas com- 
ponent, taking into account radiative cooling, heating by an 
ultraviolet (UV) meta-galactic background to emulate the 
reionization era (Haardt & Madau 1996), and a subresolu- 
tion scheme to treat SF, SN feedback and galactic winds. 
In particular, the set of parameters of the phenomenological 
wind model was chosen to obtain a escape velocity of 480 
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Figure 3. LFs for the simulated hydrodynamical galaxies at z = in comparison with the observations of Cole et al. (2001) and Huang 
et al. (2003) for the if-band (left-hand panel) and Norberg et al. (2002) for the 6j-band (right-hand panel). The simulated LFs are 
renormalized to match observations at M„ — 51og/i (indicated by an arrow; see text). Dashed and solid lines show the resulting LFs for 
the dust-attenuated and non-attenuated cases respectively. The approximate resolution limit of the simulation in each band is indicated 
by the vertical dotted lines. 



km s _1 . For the SF recipe we assume that each star particle 
of the simulation is a single stellar population (SSP) where 
the relative number of stars with different masses is obtained 
by means of the initial mass function (IMF), which in this 
work is that of Salpeter (1955) normalized in the 0.1-100 
Mq mass range. 

The metal content of gas the and stellar components 
is followed using the model presented in Tornatore et al. 
(2004; see also Tornatore et al. 2007). This scheme accounts 
for stellar evolution assessing the effect of the generation 
of Type II and Type la SN events (SNII and SNIa respec- 
tively), as well as of massive stars ending up their life in 
the asymptotic giant branch phase. These events release en- 
ergy and metals to the surrounding medium according to 
the stellar lifetimes of the different populations (Maeder & 
Meynet 1989) and to the SNIa rate (Matteucci & Recchi 
2001). It is assumed that energy feedback is only provided 
by SNs (10 51 erg per event), while metals are distributed 
in all cases following the production of Fe, O, C, Si, Mg 
and S using the stellar yields found in Recchi, Matteucci 
& D'Ercole (2001). In addition, the local cooling rate for 
each gas particle is self-consistently determined from its es- 
timated [Fe/H] abundance and temperature using the tables 
given by Sutherland & Dopita (1993). At each time step, the 
network of primordial species is solved assuming that pho- 
toionization only acts on them. In this regard, we expect 
that metal photoionization will have only a minor impact at 
the resolution covered by our simulation. 

Bound systems at the present time are identified using 
the SUBFIND algorithm (e.g. Springel et al. 2001; Dolag et 
al. 2009). Firstly, haloes having more than 32 dark matter 
particles are selected using the friends-of-friends (FoF) algo- 
rithm with a linking length of 0.2 in units of the mean inter- 
particle separation. Then, a total of 20, self bound, particles 
(ignoring gas particles) is set as a lower limit to identify sub- 
structures present in each FoF group that may be elegible 
as galaxies. As the detection of substructures is difficult at 
the low mass-end, we consider as numerically resolved galax- 
ies only substructures in the simulation having more than 



32 star particles, which results in a stellar mass resolution 
limit of ~ 7 x 10 9 Mq. 



2.1 Luminosity of galaxies 

2.1.1 Stellar population synthesis 

Similar to Saro et al. (2006) , we compute the luminosities of 
our simulated galaxies using the stellar population synthe- 
sis model of Bruzual & Chariot (2003) in different spectral 
bands. As mentioned in the previous section, each star par- 
ticle is treated as a SSP with a formation time that corre- 
sponds to redshift Zi and metallicity content Zi. For every 
galaxy, this enable us to sum up the luminosity contributions 
Li of each stellar population in the following way 

L(X) = ^2 Li(X, Zi, Zi)e~ Tx (1) 

i 

where A is the wavelength of radiation, e _TA is the extinction 
factor due to the presence of dust (see e.g. Bruzual & Chariot 
2003), and t\ is the wavelength-dependent optical depth of 
the obscuring medium (see next section). 

To compute the luminosities we use only stars inside 
the galactocentric distance that contains 83% of the galaxy's 
baryonic mass, the so-called optical radius (e.g. Nuza et al. 
2007). We have checked that adopting this condition has a 
negligible effect on the final absolute magnitudes. 



2.1.2 Attenuation by dust 

In order to asses the effect of dust in the luminosity of galax- 
ies we apply a similar procedure to Saro et al. (2009) to our 
hydrodynamical galaxies. The extinction model assumed is 
the one of Chariot & Fall (2000). These authors considered 
a scenario where the radiation of newly born stars (t < 10 7 
yr) is attenuated by the presence of dust in their birth cloud. 
Older stellar populations are able to migrate out of the 
molecular clouds and, as a consequence, are mainly affected 
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Figure 4. Halo virial velocity as a function of halo maximum cir- 
cular velocity in the simulation. The dot-dashed line corresponds 
to the best linear fit. 

by the ambient ISM. Therefore, we assume the following 
age-dependent extinction curve 

f ( TT ^y - 7 t<10 7 yr 

7-, j V 55U nm / — J 

— = { (2) 

TV \ M(^)-°' 7 *>10V 

where ry is the optical depth in the V-band, A is the wave- 
length of radiation, and fi is the fraction of the effective 
optical depth due to the ISM. In this paper we take fi ran- 
domly from a Gaussian distribution centered in 0.3 with a 
width of 0.2 and truncated at 0.1 and 1 in a similar way as 
done in previous works (e.g. De Lucia & Blaizot 2007). 

Following Guiderdoni & Rocca-Volmerage (1987) we as- 
sume that tv is proportional to the column density of cold 
gas in the line of sight having mean metallicity Z g . In par- 
ticular, for the wavelengths of interest here, we can write 

TV= (^) ( 2 .lxlo" cm-0' (3) 

where A*h is the hydrogen column density and the mean gas 
metallicity exponent s is equal to 1.6 (valid for A > 200 nm; 
see Guiderdoni & Rocca-Volmerage for details). 

To estimate the spatial dependence of -ry for our sim- 
ulated galaxies we project their gas content onto a plane 
(say xy) with the aim at averaging the cold gas distribu- 
tion as a function of projected distance to the centre. Using 
the obtained tv (r) radial profile we compute the differential 
dust-attenuated luminosity contributed by every star par- 
ticle in a given point (x, y) using its z-position as a proxy 
of its obscuration level, i.e. proportional to the amount of 
intervening cold gas between them and the galaxy outskirts. 
As a consequence, star particles behind the parent galaxy 
receive the total attenuation corresponding to its radial po- 
sition ry(r), while those ahead are almost not attenuated at 
all. The luminosity of stellar populations with intermediate 
z-positions is accordingly affected. 

As an example, in the top panel of Fig. [2] we show the 
visible optical depth averaged over projected radial bins in 
one of the most massive galaxies present in the simulation. 
As stated above, this quantity depend both on metallicity 
and column density of cold gas present in the system, hence, 
the visible optical depth tends to be higher in the inner re- 



gions of our hydrodynamical galaxies. This is due to the fact 
that, in the simulation, the galactic central regions tend to 
have higher column densities with highly enriched gas. In the 
bottom panel, we plot the fraction of luminosity attenuation 
for star particles in both bj- and A'-bands as a function of 
central distance (open and filled circles respectively). As ex- 
pected, it can be seen that the attenuation correlates with 
the optical depth profile and that it is stronger for bluer 
bands. 

Finally, we integrate over all stellar populations in order 
to estimate a global corrected magnitude for each galaxy in 
the volume (see Eq.[T|. In the rest of the paper we will always 
use dust- attenuated quantities unless explicitly stated. 



3 RESULTS 

3.1 Luminosity functions 

In Fig. Owe plot the K- and 6j-band LFs both for the sim- 
ulated galaxies at z — and for the observed local ones 
determined by Cole et al. (2001), Huang et al. (2003) and 
Norberg et al. (2002) respectively. In order to match obser- 
vations at M* — 5 log h we renormalized our model LFs in 
each case. In relation to this, we note that the number den- 
sity of L, galaxies in our simulation gives a value roughly 
two times smaller than what is observed in each band. This 
fact is also consistent with the number density we obtain 
for simulated galaxies having a typical stellar mass value 
when comparing with the local stellar mass function (e.g. Li 
& White 2009). This can be due to several reasons. Firstly, 
the background cosmology gives a baryon fraction which is 
smaller than what latest observations suggest. Secondly, the 
underlying numerical resolution adopted here (due to the 
large volume captured) can affect the assembly of small mass 
systems, thus preventing them of being the building blocks 
of more massive galaxies. These two facts combined could 
give rise to less galaxies of a given stellar mass. However, it 
is interesting to note that this is in line with the findings by 
Saro et al. (2006), where a similar deficit for galaxies within 
clusters was found, i.e. galaxy number above a given mag- 
nitude threshold tends to be smaller than observations by a 
factor of a few. 

To show the effect of our dust implementation we plot 
both the attenuated (dashed line) and non-attenuated (solid 
line) LFs. As expected, the dust extinction is more impor- 
tant in the 6j-band mainly affecting in this case to the bright- 
est (and gas-rich) systems. As can be seen in this figure, 
the simulations are not able to reproduce the bright-end of 
the LFs given by the observations. This is due to gas cool- 
ing excess onto the dark matter potential wells that leads 
to an overproduction of stars in the most massive objects. 
As mentioned in the introduction, this behavior was first 
pointed out by SAMs indicating the need of a mechanism 
capable of supressing SF as a function of time and halo mass. 
Nowadays, there is a growing body of evidence, both from 
theory and observations, that AGN feedback can play an im- 
portant role in this respect (e.g. Bower et al. 2006; Croton 
et al. 2006; Rafferty et al. 2006; Malbon et al. 2007; Kha- 
latyan et al. 2008; Lagos, Cora & Padilla 2008; Cattaneo et 
al. 2009), although the subject is far from being completely 
understood. 
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Figure 5. i\-band field LF at z = for different SF suppression models (left-hand panel: long-dashed, dotted and solid lines which 
correspond to models 2, 5 and 8 of Tabic [l] right-hand panel: three-dot-dashed, solid and dot-dashed lines which correspond to models 
7, 8 and 9 of Table [l}. Solid lines in both panels correspond to the best-fitting model. The simulated LFs are renormalized to match 
observations at M* — 5 log h (indicated by an arrow; see text). As a comparison we plot as a dashed line the dust-attenuated LF without 
SF suppression. The approximate resolution limit of the simulation in each band is indicated by the vertical dotted lines. 
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Figure 6. Idem as Fig.[5]for the fej-band field LF. 



In the same spirit of early SAMs (see below), we approx- 
imate the SF suppression in massive haloes ignoring a given 
fraction of the formed stars, using the stellar population ages 
and the maximum circular velocities of the parent galaxies 
present in the z — selected subhaloes of SUBFIND as a 
proxy. The maximum circular velocities serve as a measure 
of subhalo mass and, hence, they can be used to estimate 
the virial velocity of the systems. The relation between the 
maximum circular velocity and the virial velocity for haloes 
in the simulation can be seen in Fig. [4] For each resolved 
galaxy, we apply a selective cutoff to the stellar populations 
that were born at different epochs to mimic the required 
SF suppression and test for its effect on the LFs and other 
observables. For the sake of simplicity, we assume a linear 
relationship between the subhalo maximum circular veloc- 
ity Vrnax an d the redshift from which we do not take into 
account the newly-formed stars z C ut> as follows 

Vrnax = Vo + ZcutVl, (4) 

where Vo and Vi are free parameters. Although it is a very 



simplified parametrization aiming at reconciling the bright- 
end's LF problem, it should allow to explore the behavior 
of SF suppression in a simple way. Given a (Vo, Vr) pair, 
the net result of this recipe is to only affect systems above 
the velocity limit given by the zero point of the relation Vo, 
regulating the strength of the suppression with its slope Vi. 
The different models explored for this suppression scheme 
can be seen in Table [T] 

This procedure is similar to what was done in earlier 
semi-analytic work, where switching off gas cooling in haloes 
over a certain mass threshold was a common procedure to 
avoid excessive star formation (e.g. Kauffmann et al. 1999; 
Hatton et al. 2003; Cora 2006; Cattaneo et al. 2006). Nev- 
ertheless, it is fair to note that in our simple scheme, it is 
not possible to ignore the non-trivial interaction between 
the non-selected stars in a given halo and the surrounding 
medium for z < z cut , which will mainly affect the column 
density estimates for the suppressed systems. For instance, 
the effect of SN feedback due to the ignored stars is able 
to heat up the surrounding gas lowering the cold gas frac- 
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Figure 7. Luminosity distribution for galaxies in groups and clusters for the SDSS r- (left-hand panel) and z-bands (right-hand panel) 
in the simulation. The best-fitting to data in the RASS-SDSS galaxy cluster survey given by Popesso et al. (2006) is also shown as a 
smooth solid line. The dashed lines correspond to simulated cluster galaxies without SF quenching in massive haloes, while the solid 
ones correspond to the best-fitting SF suppression case. The simulated result is forced to match observations at M rjZ = —20 (indicated 
by an arrow; see text). The approximate resolution limit of the simulation in each band is indicated by the vertical dotted lines. 




Figure 8. Luminosity-halo mass relation for BCGs in the simu- 
lation showing the X-band luminosity as a function of their virial 
mass in comparison with the observational trend (dot-dashed line) 
given by Brough et al. (2008). Open circles stand for the SF non- 
suppression case, while filled ones represent our best-fitting SF 
quenching scenario. 

tion. In the same way, the presence of these stars in the 
simulation will also pollute the gas phase with more chemi- 
cal elements, accordingly affecting the dust properties of the 
medium. However, within the simple suppression framework 
adopted here, a more detailed model for dust extinction tak- 
ing into account these facts is beyond the scope of this paper. 

We show the outcome of applying this scheme in Figs. 
[5] and [S] for different parameter choices, where the K- and 
fej-band LFs can be seen. In the left-hand panels, three dif- 
ferent models are shown maintaining Vi fixed at 200 km 
s _1 , whereas Vb adopt the values 100, 150 and 200 km s _1 
(models 2, 5 and 8 of Table [T] long-dashed, dotted and solid 
lines respectively). As it can be seen, varying the zero point 
of the relation produces a horizontal shift in the resulting 
LF keeping its shape almost unchanged. Reducing Vb rises 
the number of suppressed galaxies thus limiting the number 



of systems at a given luminosity. In the right-hand panels, 
we also show three different models but fixing Vb at 200 km 
s" 1 and letting Vi adopt the values 150, 200 and 250 km 
s -1 (models 7, 8 and 9 of Table [T] three-dot-dashed, solid 
and dot-dashed lines respectively). In this case, all model 
LFs coincide at ~ M, — 5 log h showing appreciable differ- 
ences only in the bright-end. This reflects the strength of SF 
suppression recipe in more massive haloes. 

For the explored models, the best agreement with ob- 
servations is obtained for (Vo,Vi) = (200,200) km s _1 
(model 8, see Table[T]). This maximum circular velocity value 
roughly corresponds to Kdr ~ 150 km s _1 as can be seen 
from Fig. [4] In this model, systems with Vw r < 150 km 
s _1 are not affected by the SF quenching and those having 
Vvir > 300 km s _1 are suppressed since z > 1, in agree- 
ment with previous SAM findings which take into account 
physically motivated AGN feedback models (e.g. Croton et 
al. 2006). In what follows, we will refer to this case as the 
best-fitting model. Even though in this case there is a much 
better agreement with data in the bright-end compared to 
the non-suppression case, some improvement is still needed 
in the faint-end with the simulation results typically overpre- 
dicting the number density of galaxies per magnitude bin in 
relation to observations. As mentioned in the introduction, 
it is believed that stellar feedback should be more efficient in 
smaller galaxies. Nevertheless, higher resolution simulations 
are needed in order to study this effect. 

Given the size of the simulated volume it is possible 
to study the distribution of luminosities in the groups and 
galaxy clusters formed. However, we do not pretend here to 
make a rigorous comparison with observations since we will 
address this issue in future work using high-resolution runs. 
Instead, we just select all galaxies belonging to groups with 
virial mass higher than 10 13 hT 1 Mq as our galaxy sample. 
In Fig. [7] we show the luminosity distribution of this sample 
for the SDSS r- and z-bands together with the fit to the data 
given by Popesso et al. (2006) for the RASS-SDSS galaxy 
cluster survey (smooth solid line). To take into account the 
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Figure 9. Idem as Fig. [T] but for the galaxy distribution colour-coded from blue to red by their B — V colours (0.3 < B — V < 1) in 
the best-fitting SF suppression case. Hydrodynamical galaxies with redder colours tend to inhabit high density regions in contrast with 
bluer systems populating the field. 



(arbitrary) adopted normalization in observations our clus- 
ter LFs are forced to match the data at M r>z = —20. As can 
be seen, the best-fitting model (solid line) does a better work 
in describing the observed profiles than the dust-attenuated 
only luminosity distribution (dashed line). This is also true if 
we study the BCG luminosity as a function of halo mass. Fig. 
[H]shows the if -band luminosity-halo mass relation for BCGs 
residing in haloes with M v i r > 10 13 ft -1 Mq in both men- 
tioned cases (open and filled circles respectively), compared 
to the observational trend given by Brough et al. (2008) 
(dot-dashed line). It can be seen that the suppression sce- 
nario displays a better agreement with data, although the 
suppression seems to be somewhat stronger than needed for 
the most massive systems. 

3.2 Galaxy colours 

In Fig. [5] we show the distribution of galaxies in our simu- 
lated local volume projected onto the sky. The galaxies are 
colour coded using their B — V colours in the best-fitting sup- 
pression case. From this figure it can be stated qualitatively 
that redder hydrodynamical galaxies tend to inhabit high 
density environments, in contrast with galaxies populating 
the field, generally bluer. A similar trend is obtained for the 
plain simulation output, although in that case there is an 
excess of blue systems in cluster environments. 

The effect of SF suppression in the colours of galax- 
ies is easily seen in Fig. 1101 There we show the B — V 
colour index of numerically resolved galaxies against their 
stellar masses for the non-suppression case (filled circles), 
and what results for the best-fitting suppression procedure 
presented above (contour lines). Interestingly, the hydrody- 
namical galaxy population of the non-suppresed simulation 
shows a clear colour bi-modality for a given stellar mass, 
where two main branches can be easily seen. Those systems 
with colour index satisfying B — V > 0.8 populate the so- 
called red sequence, while the rest populate what we call 
the blue cloud. This last branch contains a blue galaxy pop- 



ulation with M, > 2 x 10 11 h~ x M responsible for the 
mismatch in the bright-end of the LFs shown in Fig. [3] In 
particular, most of the systems having V v i T > 300 km s _1 
(arrow in the right) correspond to overluminous BCGs which 
have acquired around 50% of their stellar mass at late epochs 
due to the overcooling effect. 

This is at variance with the SF suppression scenario 
(contour lines) in which most of the massive blue galaxies 
became redder as a result of their older stellar populations. 
It can be seen that many of the systems with V v i T > 150 
km s _1 (arrow in the left) have been suppressed in such a 
way that their colours turn redder than what observations 
suggest. Henceforth, it is clear from the contours that the 
simple procedure adopted to suppress SF is too efficient for 
some galaxies, while not efficient enough for the most mas- 
sive ones. This results in a distribution that cannot fully ac- 
count for the galaxy colour bi-modality in the stellar mass 
range M, ~ 2 x 10 10 - 2 x 10 11 h' 1 Mq (e.g. Baldry et 
al. 2006), reflecting the fact that physically motivated feed- 
back schemes must be considered for massive haloes in future 
work. 

As a final remark, it is worth noting that a moderately 
massive red population having M* ~ 5 x 10 9 h -1 Mq nat- 
urally appears in the simulation. The diamond in Fig. [10] 
indicates the approximate position of this population in the 
diagram. These systems correspond to galaxies that have 
assembled most of their stellar mass at early epochs, being 
satellites which gas content have been stripped off when ac- 
creted by central objects. This mechanism removes gas from 
galaxies limiting new SF and, as a consequence, stellar pop- 
ulations become older and redder. 



3.3 Clustering properties 

In this section, we explore the clustering properties of galax- 
ies in the simulation, both for the non-suppresed and best- 
fitting SF suppression cases. The left-hand panel of Fig. 
1111 shows the real-space two-point correlation function for 
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m, [h- 1 M ] 

Figure 10. B — V galaxy colours as a function of their stellar 
mass for the numerically resolved galaxy population at z = 0. 
Filled circles stand for the case without SF quenching, while con- 
tour lines contain 55%, 85% and 95% of the galaxy population 
in the best-fitting SF suppression scenario. Arrows indicate lower 
bounds for the virial velocity of systems with a given stellar mass: 
V v i r > 150 and 300 km s _1 respectively. The diamond indicates 
the approximate position of the satellite galaxy population in the 
diagram. 



ulation makes the number density of red systems low, thus 
rising the value of the correlation length. 

On the other hand, Fig. 1121 shows the outcomes when 
considering the best-fitting suppression scheme. From its left- 
hand panel, it can be seen that the correlation function for 
~ Z/» galaxies (circles) exhibits essentially the same behav- 
ior as before (compare with Fig. Hip , while displaying ap- 
preciable differences only for brighter galaxies (squares and 
diamonds). For objects brighter than Mb i — 5 log ft = —21 
the result is affected by low number statistics leading to a 
noisier correlation function that is somewhat off from the 
observational expectation. This is partially overcome when 
including more galaxies in the M&, — 5 log h > —20 case. Ac- 
tive galaxies present an excellent agreement with data for all 
studied scales (right-hand panel of Fig. I12|l , and the discrep- 
ancy with observations found for passive red systems found 
in Fig. [TJJ is alleviated. This is somewhat expected since, 
as can be seen from the contour lines of Fig. 1101 after the 
SF quenching, many of the galaxies have migrated to the 
red sequence rising the number density of passive systems. 
Despite of this improvement, it is important to remember 
that the colour-stellar mass relation is not well reproduced 
in this case due to the simplistic approach adopted to redden 
galaxies more massive than Vvir > 150 km s _1 . 



galaxies in the plain simulation with different fy-band ab- 
solute magnitudes (filled symbols). Circles are related to 
galaxies with —20 < Mb, — 51og/i < —19, while squares 
and diamonds show the resulting correlations for all galax- 
ies brighter than —20 and —21, respectively. We also plot as 
dashed and three-dot-dashed lines the observational trends 
found by Norberg et al. (2001) for ~ L* and ~ 4L t galaxies 
(corresponding to Mb, — 51og/i ~ —19.5 and —21.5 respec- 
tively). In the first place, it can be seen that brighter galaxies 
cluster strongly than fainter ones, as expected. The agree- 
ment with observations is fairly good taking into account the 
fact that normalization is a natural outcome of the simula- 
tion. In particular, for ~ L* galaxies (circles) the clustering 
pattern shows a power-law behavior similar to that observed, 
although the agreement is worse for higher scales (> 10 h^ 1 
Mpc). Likewise, for systems with Mb- — 51og/i < —21 (dia- 
monds) the agreement with observations is encouraging. 

It is also interesting to investigate the clustering of 
galaxies as a function of their SF activity parametrized us- 
ing galaxy colours. At z ~ 0, passive systems tend to pop- 
ulate high density environments, having redder colours, in 
contrast to active systems which are mainly located in the 
field (e.g. Balogh et al. 2004; Baldry et al. 2004, 2006). We 
define as passive systems all galaxies having B — V > 0.8, 
while active systems satisfy B — V < 0.8. The right-hand 
panel of Fig. [TTJ shows the real-space two-point correlation 
function for passive red galaxies (diamonds) and bluer ones 
(circles) in comparison with the observed correlations (valid 
for scales between ~ 0.1-20 Mpc) given by Madgwick et 
al. (2003). It is easy to see that both populations are clearly 
separated, with redder systems clustering more strongly, and 
displaying a similar behavior than observations. For passive 
systems, there is a substantial excess of power and a clear 
departure from a power-law. The power excess in this case 
is due to the fact that the overcooling present in the sim- 



4 SUMMARY AND CONCLUSIONS 

In this paper, we have focused on the luminosity, colour 
and clustering properties of galaxies that arise naturally in 
a hydrodynamical cosmological simulation of structure for- 
mation. The spherical local volume of ~ 130 3 h~ 3 Mpc 3 was 
simulated using the zooming technique with more than ~ 50 
million gas and dark matter particles. In particular, our sim- 
ulation code include different physical mechanisms such as 
supernova feedback, galactic winds, an UV pervading back- 
ground since z ~ 6 and chemical enrichment of the gas and 
stellar populations as time elapses. However, we have not 
explicitly included a physical mechanism capable of quench- 
ing SF in massive haloes (like AGNs). In this way, by z — 
a total of ~ 20000 objects having more than 32 star parti- 
cles are formed in the volume. The luminosities of the hy- 
drodynamical galaxies were computed using the differential 
contribution provided by every star particle present inside 
the optical radius of the systems. As a way of estimating 
the dust-attenuated luminosity in the bands of interest, we 
implemented a simple model for dust extinction that makes 
use of the cold gas distribution in each galaxy as a tracer 
of its dust content. In order to mimic the required SF sup- 
pression in massive (i.e. VW > 150 km s _1 ) haloes we also 
applied a simple recipe to the present day output of the 
simulation to reconcile the bright-end of the simulated LFs 
with observations and to further investigate its inclusion on 
other properties of the simulated galaxies. In the following 
we summarize the main conclusions of the present work: 

• As shown earlier by semi-analytic work, it is essential 
to suppress SF in massive haloes to describe the bright-end 
of the field LF at z ~ 0. For instance, haloes with masses 
corresponding to Vvir ~ 300 km s~ x are required to quench 
their SF since z ~ 1 in the best-fitting SF suppression 
case. In this way, it is possible to avoid the formation 
of very bright blue BCGs in the simulation. However, 
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Figure 11. Real space two-point correlation function for the simulated galaxies without SF quenching (filled symbols) in different 6j-band 
magnitude bins (left-hand panel) and colour populations (right-hand panel). The error bars denote Poissonian fluctuations. Observational 
trends from Norberg et al. (2001) for ~ L* and ~ 4L* galaxies, and from Madgwick et al. (2003) for active and passive galaxies are 
shown as dashed and three-dot-dashed lines respectively. The observations are valid in the range ~ 0.1-20 h^ 1 Mpc. 




Figure 12. Idem as Fig. 1111 for galaxies in the best-fitting post-processing scheme to supress SF. 



after a small renormalization to L*, the faint-end of the 
simulated field LFs tends to overpredict the number density 
of galaxies per magnitude bin in the 6j-band for galaxies 
with luminosities lower than ~ L*. To better study this 
effect higher resolution simulations are needed. 

• The group/cluster bright-end luminosity distributions 
are better described when SF suppression is considered. 
This is due to the fact that massive cluster satellites used 
to be BCGs in the past, and also require some level of SF 
suppression. However, as a consequence of our relatively 
high-mass galaxy resolution limit we cannot confirm the 
upturn of faint cluster satellites in the cluster LF function 
as previously obtained by Saro et al. (2006, 2009) using 
high-resolution resimulations of galaxy clusters. On the 
other hand, the Jf-band luminosity-halo mass relation that 
naturally arises in the simulation also needs further SF 
quenching in massive haloes to better describe observations 
given by Brough et al. (2008). 



• The stellar masses of the hydrodynamical galaxies in 
the plain simulation display a clear B — V colour bi-modal 
behavior as suggested by observations (e.g. Baldry et 
al. 2006), although the active blue galaxies present very 
massive systems, due to overcooling, contrary to what 
data shows. When taking into account the SF suppression 
scheme, the B — V galaxy colours turn redder, displaying 
a colour- stellar mass distribution that is not able to fully 
account for observations. In this case, a bi-modal distribu- 
tion is observed for stellar masses lower than ~ 2 x 10 10 
h^ 1 Mq, while for higher masses most of the galaxies are 
red. This is a manifestation of the rough nature of the sim- 
ple parametrization adopted to quench SF in massive haloes. 

• The clustering properties of different hydrodynamical 
galaxy populations that naturally arise in the simulation 
are in qualitative agreement with observational results. The 
scale-dependent correlation function of our ~ L* objects in 
the 6j-band displays a power-law behavior up to ~ 10 h~ 1 
Mpc that is in good agreement with data given by Norberg 
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ct al. (2001). As expected, higher luminosity systems 
show a stronger clustering pattern, where galaxies having 
luminosities around ~ 4L* display an excellent agreement 
with observations. On the other hand, it is also possible 
to disentangle the galaxy correlation properties when 
the sample is splitted by galaxy colours. The correlation 
function of active blue (B — V < 0.8) systems is in rough 
agreement with the results given by Madgwick et al. (2003). 
However, passive red systems (B — V > 0.8) present an 
excess of power in comparison with observations, showing 
a clear departure from the power-law behavior. As shown 
when using the SF suppression scheme, this discrepancy is 
alleviated when more galaxies are able to populate the red 
sequence of the colour-magnitude diagram. This is due to 
the fact that, typically, the added systems tend to reside in 
less clustered haloes, thus lowering the correlation length 
value. 

In summary, we conclude that analyzing the global prop- 
erties of the galaxy population within hydrodynamical, cos- 
mological simulations, start to be a promising tool to study 
galaxy evolution. Current simulations already fairly repre- 
sent the underlying hydrodynamical effects (at least in a 
global sense) and, in general, describe the star formation 
process well enough to qualitatively reproduce observed en- 
vironmental trends. However, similar than for the widely 
used SAMs, overcooling in massive halos has to be quenched 
by additional feedback effects. It has to be seen in future sim- 
ulations, which directly include such additional processes, if 
such a quenching happens mildly enough not to destroy some 
of general trends already captured in the current generation 
of hydrodynamical simulations. 
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